Verifying the Kugo— Ojima Confinement Criterion in Landau Gauge Yang— Mills theory 



Peter Watson and Reinhard Alkofer 

Institute for Theoretical Physics, Tiibingen University, Auf der Morgenstelle I4, D-72076 Tubingen, Germany 

(February 1, 2008) 

Expanding the Landau gauge gluon and ghost two-point functions in a power series we investi- 
gate their infrared behavior. The corresponding powers are constrained through the ghost Dyson- 
Schwinger equation by exploiting multiphcative renormalizability. Without recourse to any specific 
truncation we demonstrate that the infrared powers of the gluon and ghost propagators are uniquely 
related to each other. Constraints for these powers are derived, and the resulting infrared enhance- 
ment of the ghost propagator signals that the Kugo-Ojima confinement criterion is fulfilled in 
Landau gauge Yang-Mills theory. 
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Modifications to the standard framework of quantum 
field theory are necessary to accommodate the confine- 
ment phenomenon in QCD. Here either relaxing the prin- 
ciple of locality or abandoning the positivity of the repre- 
sentation space (or both) is required. As quantum gauge 
field theories in covariant gauges demand indefinite met- 
ric spaces, and as locality is of fundamental importance, 
it is suggestive to relate confinement to the violation 
of positivity in Yang-Mills theories: No colored states 
should be present in the positive definite space of physi- 
cal states defined by some suitable condition maintaining 
physical S'-matrix unitarity. Exploiting the BRS symme- 
try of the gauge fixed action a nilpotent BRS-charge 
in an indefinite metric space V can be constructed. The 
subspace Vp = Ker Qg is defined by those states which 
are annihilated by the BRS charge Qg. Positivity is then 
proved for physical states ^ (see also refs. [§|j3|) which 
are given by the cohomology H{Qb, V) = Ker Qg/ImQg, 
the covariant space of equivalence classes of BRS-closed 
modulo BRS-exact states. 

Longitudinal and timelike gluons form BRS quartets 
together with ghosts and are thus unphysical. At the 
same time the global symmetry J° corresponding to 
gauge transformations generated by 9°'{x) = aj^x'' is 
spontaneously broken quite analogous to the displace- 
ment symmetry in QED. The identification of BRS- 
singlet states with color singlets is non-perturbatively 
possible if this global symmetry is restored The 
condition for a dynamical restoration of this symmetry is 
the Kugo-Ojima confinement criterion which is a require- 
ment on the infrared behavior of a four-point Green's 
function. In Landau gauge, a sufficient condition for this 
criterion is that the nonperturbative ghost propagator is 
more singular than a massless pole in the infrared M : 
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This mechanism is correlated to the derivation Q of 
the Oehme-Zimmermann superconvergence relations ^] 
from Ward-Takahashi identities. These superconver- 
gence relations formalize a long known contradiction be- 



tween asymptotic freedom and the positivity of the spec- 
tral density for transverse gluons in the covariant gauge. 

The violation of positivity of transverse gluons is un- 
ambiguously established by a variety of independent non- 
perturbative studies of the gluon propagator 0| . Further- 
more, recent lattice calculations investigate the Kugo- 
Ojima criterion directly Q. As infrared singularities are 
anticipated, also non-perturbative continuum methods, 
besides lattice calculations, are necessary to investigate 
this and related pictures of confinement in more detail. 

The Green's functions of a quantum field theory obey 
the Dyson-Schwinger equations. For realistic theories 
this infinite tower of integral equations cannot be solved 
directly, and truncations are necessary to obtain (numer- 
ical) solutions. The ones for QCD propagators can then 
be used to formulate a successful hadron physics phe- 
nomenology, see e.g. the recent reviews 

The results of two recent truncation schemes ^0 11 
suggest that the Kugo-Ojima confinement criterion for 
the ghost propagator and the violation of positivity for 
transverse gluon correlations are fulfilled. Obviously, the 
question arises whether these results are sensitive to the 
employed truncation schemes. 

In this letter we will exploit the general structure of the 
simplest QCD Dyson-Schwinger equation, the one for the 
ghost two-point function, to relate the infrared behav- 
ior of gluon and ghost propagators. Thereby the Kugo- 
Ojima confinement criterion will be verified by purely an- 
alytical methods. Using the rather unique properties of 
the ghost-gluon vertex in the Landau gauge and assum- 
ing the applicability of a powerlaw ansatz in the infrared 
the ghost Dyson-Schwinger equation admits a completely 
general qualitative analysis. Instead of truncating the 
system, the general forms for the unknown functions con- 
strained by various consistency arguments, notably mul- 
tiplicative renormalizability, are input into the equation 
under an asymptotic expansion in the infrared momen- 
tum scale, c.f. the method described in ref. fl^ . 

We start by providing some basic notations. The gen- 
eral form of the gluon propagator in Landau gauge is 
purely transverse and is given by 



(2) 



which defines the gluon renormalization function Z(j)^). 
The ghost renormaUzation function G{p'^) is estabhshed 
in eq. (|l|). We define the ghost-gluon vertex 

Tf-{pi,P2,P3) = -»5r^%(pi,P2,P3) (3) 

with antighost (pi), ghost {p2) and gluon (pa) momentum 
incoming, i.e. pi +P2+Pz = 0. At tree-level the dressing 
functions reduce to 



Z{p^) - G(p2) ^ 1, r^(pi,P2,P3) =Pl, 



(4) 



These renormalized dressing functions are related to the 
bare ones by the usual renormalization constants 
(gluon wave function), Z3 (ghost wave function) and Zi 
(ghost-gluon vertex). In Landau gauge one has Zi = 1 
[^,0, i.e., the ghost-gluon vertex has not to be renor- 
malized . 

The renormalized ghost dressing function G(p^) is a di- 
mensionless function of a single variable. In pure Yang- 
Mills theory there are only two scales present: the exter- 
nal scale p^ and the renormalization scale /i^ . Therefore 
G is a function of the ratio p^ / only. As the external 
scale vanishes this function can be written as an asymp- 
totic expansion. Renormalizability implies that it must 
also depend on the renormalized coupling g. This moti- 
vates the power series ansatz 



G{p'■^l'') = Y^d 2 



(5) 



where the coefficients dn and the exponents 5n are, at 
least in principle, functions of g. We will see below that 
the exponents (S„ will be independent of o. Furthermore, 
as will be shown later on, the series (|^) has a lowest 
power i5o such that as the external scale p^ vanishes the 
ghost dressing function G{p'^) is well approximated by 
a powerlaw. Multiplicative renormalizability allows to 
derive the renormalization group (RG) equation 



M|^+/3(5)^-27G)G(p2;/i') = (6) 



with /3(g) being the beta function, and 

dXogZ--^ 



7g = M- 



(7) 



is related to the ghost anomalous dimension. Inserting 
the series (||) into eq. (|^) , and exploiting that p^ is a free 
(though small) parameter one obtains 

- d,An + I3{g) f ^ + d„ log(^)^ j + rf„7G = 0. (8) 



Since the logarithm has no matching factor the expo- 
nents Sn must be independent of the renormalized cou- 
pling {d6n/dg = 0). Using (3{g) = -5(270 -I- 7a) with 
_ j_ 

7a ~ fid log Z^ ^/dfi the resulting differential equation 
for the coefficients d„ has the general solution 



dn = const, g ^ic+fA . 



(9) 



Expanding the renormalized gluon dressing function 

„2 \ 



Z{p ;fi 



one arrives in a completely analogous way at 

2(€n + Tr,) 

e„ = const, g ^■'c+tA . 



(10) 
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In the renormalization process we have to specify the 
corresponding normalization conditions. We require that 
the unrenormalized functions reduce to a constant at 
asymptotically large momenta. Exploiting additionally 
multiplicative renormalizability and using eq. (^) one ob- 
tains the formal relation 



Z3(/.2,A2)0C [Y^dn 



A^ 



(12) 



This equation is only well-defined if the corresponding se- 
ries makes sense, e.g. if it is resumable by some method. 
This is the main assumption throughout the presented 
derivation. One consistency check is provided by insert- 
ing the general solution for the coefficients d„ (^ and 
rewriting the renormalized coupling g in terms of the 
scale fi. This yields 



Z3(/i',A2) cx 



7G 



(13) 



which is precisely the correct dependence to agree with 
the definition of 7g. A completely analogous argument 
applies to the gluon renormalization constant Z3. 
The tensor structure of the ghost-gluon vertex reads 



r^(Pl,P2,P3) ^Plf^MPi) +P3t,B{pf) 



(14) 



In Landau gauge the ghost-gluon vertex becomes bare 
when the ghost momentum p2 vanishes ||l3|,0 : 



^^liPl,P2,P3) =PiM (^ - B) -p2t_iB 



P2^0 



(15) 



Thus, as p2 — > the functions A and B cannot be singu- 
lar. This is in accordance with the fact that the renor- 
malization constant Zi = 1. 

The RG equation for the vertex function can be split 
into two parts, one for each of the tensor components. 
It is only necessary to consider the equation for ^4, and 
since Zi — 1, it is identical to the RG equation for the 



coupling. Writing the vertex function as a series has the 
comphcation of several momentum scales. For our con- 
siderations it is sufhcient to simply extract the scale and 
coupling dependence 



k 



with 



(16) 
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We must renormalize the ghost-gluon vertex at the point 
P2 = since this corresponds to a zero finite renormal- 
ization from Zi = 1. Knowing the P2 — > limit of the 
full vertex does not constrain A. (When P2 = one only 
has A — B — 1.) With this given limit and the additional 
requirement that the ghost equation should provide a 
consistent solution one infers that 



k 



(18) 



is a finite positive constant. Asymptotic freedom requires 
that 7a + > 0. In the unrenormalized theory the 
vertex should be bare when the unrenormalized coupling 
vanishes because Zi = 1. Therefore Kk < with the 
highest value kq = 0. 

The ultraviolet region of the ghost Dyson-Schwingcr 
equation corresponds to a small (albeit non-vanishing) 
ghost momentum p2 and all other momenta large. For 
the following discussions it is helpful to write 



k.7n 



bk{p\)>^k 



(19) 



where the p„i > and po = 0. 

The full Euclidean ghost Dyson-Schwinger equation 
can be written after having performed the two trivial 
angular integrals as 



G-\x)^Z:i 

/ dyZ{y)- 
IGtt^ Jo Jo 



(20) 



d0sm^0^^A{z,x,y) , 



where 9 is the angle between the two internal loop mo- 
menta corresponding to y = and z = {p — q)^ with 
X = p^. Ca = Nc is the resulting color factor. When 
using the expansions (||) and ( [lo| ) the angular integra- 
tion splits the radial integral into two parts, the lower 
part Igh (0 x) and the upper part Jgh {x — > A^). The 
expression derived from the upper limit will be some com- 
bination of factors with the ratio x/A^ to some power. 
The effect of an exact angular integration is to generate 
a hypergeometric function whose argument is x/y and 
which can be expanded in a convergent series with posi- 
tive integer powers of x/y starting with unity. Thus, the 
^-dependence of the vertex function is unimportant. 



The most important aspect is that the ghost-gluon ver- 
tex reduces to its tree-level form as the in-ghost momen- 
tum vanishes. In the ultraviolet part of the integral this 
is precisely the case: the in-ghost momentum is the in- 
frared external scale x. Therefore the vertex function can 
be written in the form of eq. (|l9|). 

Using the series for G and Z, it is possible to write the 
ghost equation with the most general vertex function as 




= Z. 



(21) 



CaO^ XT' ^ I 

72 1^ Sj^fc [ -2 

J,k 



Ighiej,Sk) + Jgh{ej,Sk) 



where Igh and Jgh represent the lower and upper regions 
of the integral, respectively. The integrals Igh must con- 
tain factors characterized by the Kk but there are no re- 
strictions on its form except that the vertex must have 
positive powers of the coupling. As the integral Igh is 
very difficult to compute we will not be able to proceed 
by studying its general properties. 

It is frequently asserted that only the infrared inte- 
grals of Dyson-Schwinger equations can contribute to the 
renormalized propagator functions on the left hand side 
because this is the only fully renormalizable possibility. 
Indeed, one can show that upper integral Jgh will not 
contribute directly to the propagator functions |l^ . It is 
the renormalization constant Z3 that provides the key to 
the infrared behavior of the propagators, and only the ul- 
traviolet integral can yield the ultraviolet divergent part 
of the ghost equation. First we note that Jgh{ej,Sk) can 
be expanded as 



•^gh — Jj,k,l,m,q [ ^ 



2 \ i^j+Sk-Pm-q 
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, (22) 



where pm > represents the kinematical content of the 
vertex function, and q = 0,1, . . . denotes the order of the 
angular integration series. As Z3 has to be independent 
of x one deduces from eqs. (^,^ that only po = q = 
will contribute. All other terms involve x/A^ to a positive 
power and subsequently vanish: 



^3(m',A2) 

167r2 



3-k,l 



ejdkJj.kA.o.o^i 



^2 \ <ij+Sk+K.l 



(23) 



Using that expression (|18|) is a constant, employing the 
expansion ( [T^ ) and multiplying through with the left 
hand side yields 



^ejdidkJo.o,Q[£],Ok) I — 
V P 
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i,i,k 



const. 



(24) 



Since eo and Sq are the unique lowest powers, and the 
term with i = j = k — has nothing to compete with 
it must be responsible for the non-vanishing part of the 
left hand side of this equation. Therefore it must be 
independent of A^//i^, and wc obtain 



eo = -2(5o 



(25) 



i.e. we have shown that the leading infrared powers of 
the gluon and the ghost propagator are uniquely related. 
Note that this result is valid even if corresponding series 
are not convergent. We have simply matched the coeffi- 
cients of an asymptotic series. Therefore, our underlying 
assumption is that the employed series have to exist in 
the sense of being, at least, asymptotic ones. 

Eq. (Q) constrains the values of the related leading 
powers. To do this, one needs to explicitly evaluate the 
integrals Jj^k,o,o,o which are readily calculated from the 
full ghost equation with a bare vertex stripping off the 
terms that depend on the external scale x. Concentrating 
on the leading powers we obtain: 



1 



> 0. 



(26) 



To yield a positive ghost renormalization function {i.e. to 
keep the negative-norm property of the ghost states) the 
corresponding leading order coefficients eo and do niust 
be positive. As the coupling is positive this yields: 



<5n < 



(27) 



This relation implies that the Kugo-Ojima confinement 
criterion^ see eq. (|^), is fulfilled. 

A lower bound on 5q can be derived from the lower 
bound of the integral in the ghost equation. As this 
bound is considerably weaker as the one following from 
a generalized masslessness condition (see eq. (192) in ref. 
1^) (5o > — 1, we will not present its derivation here. This 
implies 2 > eo > 0: The gluon propagator is suppressed 
in the infrared as compared to its tree-level form. 

In summary, we have investigated the general structure 
of the ghost Dyson-Schwinger equation in Landau gauge. 
A general series form for the propagator and vertex func- 
tions has been employed. No specialized ansatz for the 
ghost-gluon vertex has been used. Care has been taken 
not to use approximations that destroy the generality of 
the results. Multiplicative renormalizability, and espe- 
cially the non-renormalization of the ghost-gluon vertex, 
have placed the constraints to relate uniquely the leading 
infrared behavior of the gluon and the ghost propagator. 

We have found an infrared enhanced ghost propagator 
thereby verifying the Kugo-Ojima confinement criterion 
|0,D directly. Due to the relation ( p5| ) this corresponds 
to an infrared suppressed gluon propagator. For a non- 
integer leading infrared power this implies that positiv- 
ity is violated and the Oehme-Zimmermann superconver- 
gence relation is fulfilled non-perturbatively. Consid- 
ering only the ghost equation as in this letter an infrared 



constant gluon propagator cannot be excluded. To inves- 
tigate this question one has to consider additionally the 
gluon Dyson-Schwinger equation |]l5| . 
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